or 



> 

in 



o 
oo 



- '-I 

>< 



Quantum Computation in a radio-single mode cavity: the dis- 
sipative Jaynes and Cummings Model 



OO ! P.Ferrante 

o 



Dipartimento di Ricerche Energetiche ed Ambientali, Universitd degli Studi di Palermo 

(N 

S,Nicolosi 
Vh ■ 

Dipartimento di Ricerche Energetiche ed Ambientali, Universitd degli Studi di Palermo 

■^T CNR-LATO hsr Giglio, CefaliiB nicolosiQdream.unipa.it 



G.Scaccianoce 

Dipartimento di Ricerche Energetiche ed Ambientali, Universitd degli Studi di Palermo 

G.Rizzo 

Dipartimento di Ricerche Energetiche ed Ambientali, Universitd degli Studi di Palermo 

a 

O^ (Received 2007) 



Abstract. Our aim is to give an account, trough an analysis of a number of papers by 
F.Petruccione, H.J.Charmicael, J.C. Raimond and their contemporaries, of the specific answer 
that we gave to the problem of Open Quantum Systems dynamical evolution and how this idea 
evolves and develops in physical research and in the scientific debate of the following decades. 

f/-\ , Permanent solution should not been accepted from physical research, but analysis of the real 

work of scientists, of the difficulties they face and the ever changing solutions they offer is, we 
believe, part of our understanding of science and an indispensable basis for further methodological 
inquiries. 

In this paper we have chosen to analyze a dissipative Jaynes-Cummings model assuming the 
common electrodynamics free field for the bipartite system and an another independent bath for 
the cavity, so taking into account loosing of energy because of the imperfect mirrors. The Nud 
theorem application leads to predict new cooperative effects between the atom and the cavity 
mode as the creation of conditional transient entanglement, tending to become stationary as the 

r » "] ■ coupling constant take a well defined value. 

1. Introduction 

The idea of implementing quantum information devices based on the use of 
single atoms or molecules has gone progressively growing up in the course of the 
last few years. The reason can be envied in the high contemporary evolution of 
theoretical dynamical models along with the ability reached by experimentalists in 
manipulating quantum objects first considered theoretician's tools pQ. The state 
of art is that, although it has been possible to obtain entanglement conditions 
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between elementary systems in different physical scenarios [I], the temporal per- 
sistence of quantum coherence is an open problem. In this moment, it is therefore 
necessary to put attention on the theoretical aspects of decoherence in order to pre- 
dict the experimental conditions under which it is maximally reduced [3J. Almost 
Decoherence Free Substates seem to have the characteristics of eligibility needed to 
implement quantum computation. Their generation and temporal persistence can 
be theoretically predicted in some high symmetrical models [U \7\ EJ QUI E] • 
The application of the early introduced formal solution of the Markovian Master 
Equation (Nud theorem) [3j has supplied, in the few analyzed systems, the pre- 
diction of a conditional building up of entanglement. The result is obtained under 
symmetrical condition corresponding to specific locations of involved subsystems 
(atoms, molecules). Unfortunately, despite the positive results in isolating single 
quantum objects [I], the difficulties connected with the location of more than one 
atoms in fixed arbitrary points is an open task. 

In general, practical quantum gates implementation needs to meet stringent 
criteria in order to operate successfully [5]. First of all, the qubits has to be suf- 
ficiently isolated in order to manipulate them in a controlled environment. They 
need to be initialized precisely. The effective interactions among qubits should be 
carefully tuned and a set of quantum operations should be made possible in order 
to perform any other required quantum gates. Finally, the system must be scalable 
to more than a few qubits |12j . 

These considerations have encouraged us to speculate about the possibility to 
use hybrid entanglement conditions in order to implement quantum protocols. A 
single mode cavity and a two level atom may be considered a hybrid two qubits 
system. The isolation of a single atom for enough long time inside a cavity is today 
possible. Moreover, the Jaynes and Cummings Model is one of the few exactly 
soluble models in quantum optics. It predicts many interesting non classical states 
experimentally testable in laboratory [1\. In realistic situation, however, one per- 
forms experiments in cavity with finite Q and in presence of atomic spontaneous 
emission. So it become of fundamental importance to know how the predictions 
of this model are affected by the unavoidable presence of loss mechanisms. This 
problem is currently very extensively studied, the main approximation being the 
assumption of two different reservoirs, one for the atom and one for the cavity 
mode respectively [2]. The further calculation are simpler under this approxima- 
tion because the two bath don't introduce coherences between the atom and the 
cavity: the master equation is in Lindblad form and it is easier to solve than the 
one we derive; the solution being the total destruction of coherences because of 
the two different and nonspeaking channels of dissipation. Moreover, my equation 
contains the simple one as particular case. Actually, here we show to be able to 
find the exact solution of a dissipative J-C model assuming a common reservoir 
for the bipartite system, which, on the ground of the above consideration, appears 
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to be a more realistic hypothesis. This leads to the prediction of new cooperative 
effects, induced by the zero-point fluctuations of environment, between the atom 
and the cavity mode as the creation of conditional transient entanglement, tend- 
ing to become stationary as the strengths of the coupling with the reservoir take a 
well defined value. Finally, in order to be maximally realistic, Prof. J.M.Reimond 
suggested me to consider also the loose of energy due to the imperfect reflection 
on the mirrors. This correction in the microscopic model does not introduce com- 
plication in the solution of the relative master equation because the new bath is 
independent from the first one and, indeed, easily treating from a theoretical point 
of view (not induced coherences). In presence of the second channel of decoherence 
the building up of entanglement exists during the transient period in which the 
atom is confined inside the cavity. The long time solution (the order of magnitude 
of the time involved is given in the next sections) shows that the introduction 
of the second bath makes disappear coherence (Rabi oscillation) among the two 
subsystem involved: the atom and the cavity mode. Despite this fact, the time 
involved in decoherence process can be made much longer than those necessary 
to implement a quantum protocol as deducible from the theoretical analysis here 
developed if interfaced with the experimental measures performed by Aroches' 
group. The measures appear well fitted by the theoretical model here proposed. 
The standard models (two different baths) are able to reproduce only the top of 
the curve (dissipation [4j). Instead mine is able to reproduce also the lower part 
(cooperation [1]). 

The paper is structured as follow: in section II we report the principal step 
and approximation leading to the microscopic derivation of the Markovian Master 
Equation (MME) and we solve it when T = 0, showing also the full equivalence 
between MME and Piecewise Deterministic Processes (PDP) [3]. In section III 
we apply the NuD theorem to derive the solution to the dissipative J-C model. 
In section IV we try to justify the obtained dynamical behaviour in terms of 
continuous measurement theory. 

2. The Markovian Master Equation 

It is well know that under the Rotating Wave and the Born-Markov approxima- 
tions the master equation describing the reduced dynamical behavior of a generic 
quantum system linearly coupled to an environment can be put in the form [4] [13] 

p s (t) = -i[H s + H LS , p s (t)\ + D(p s (t)), (1) 

where H$ is the Hamiltonian describing the free evolution of the isolated system, 

w a, 13 

-\{ASv)MuUs(t)}), (2) 
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Hls = EE^M4MA^), (3) 

S a ^) = ±.{T a ^)-Tl a (u)) (4) 

and 

7a,/3(w) = r ai/3 (u0 + rj >a (w), (5) 

r a ^(a;) being the one-sided Fourier transforms of the reservoir correlation func- 
tions. Finally we recall that the operators A a (u>) and A< a (u>), we are going to define 
and whose properties we are going to explore, act only in the Hilbert space of the 
system. 

Eq. (TjQ) has been derived under the hypothesis that the interaction Hamiltonian 
between the system and the reservoir, in the Schrddinger picture, is given by [13] 

fl> = E^ a ®S a , (6) 

a 

that is the most general form of the interaction. 

In the above expression A a = (A a y and B a = (B a y are operators acting 
respectively on the Hilbert space of the system and of the reservoir. The eq. 
dU) can be written in a slightly different form if one decomposes the interaction 
Hamiltonian into eigenoperators of the system and reservoir free Hamiltonian. 

DEFINITION 1. Supposing the spectrum of H$ and Hb to be discrete (general- 
ization to the continuous case is trivial) let us denote the eigenvalue of H $ {Hb) by 
s (77) and the projection operator onto the eigenspace belonging to the eigenvalue 
e (rj) by n(e) (11(77)). Then we can define the operators: 

A a (u;)= E n(e)A Q n( £ '), (7) 

s —e=w 

B a (u)^ E n( V )B a II( V '). (8) 

/ 

T] —77=0) 

From the above definition we immediately deduce the following relations 

[H S , A a {u)\ = -uAaiu), [H B , B a {u>)] = -uB a {u), (9) 

[ff s ,4( W )] = + W 4H and [H B ,Bi{uj)]=+u J Bi(uj). (10) 

An immediate consequence is that the operators A^uj > 0) e A a (uj > 0) raise 
and lower the energy of the system S by the amount Tiuj respectively and that the 
corresponding interaction picture operators take the form 
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e iHst A a {uj)e- iHst = e~ iu!t A a {uj), e iHBt B a {uj)e- iHBt = e~ iult B a (u;), (11) 

e iHst Ai(Lo)e- iHst = e +i ^Ai(io) and e iHBt Bl (oo)B- iHBt = e +iu)t B^uj) . (12) 

Finally we note that 

Al{uj) = A a {-uj) and B^u) = B a (-u). (13) 

Summing eq. (|13p over all energy differences and employing the completeness re- 
lation we get 

£4,(a;) = £4x(-w) = 4, and J2 B l ( w ) = E B «(-w) = B a (14) 

LU LU LU LU 

The above positions enable us to cast the interaction Hamiltonian into the following 
form 

Hi= Y. A a {uj) ® B a {J) = £ Al{u)®Bi{J). (15) 

The reason for introducing the eigenoperator decomposition, by virtue of which 
the interaction Hamiltonian in the interaction picture can now be written as 

Hi(t)= J2 e-^+^A^^B^u;'), (16) 

is that exploiting the rotating wave approximation, whose microscopic effect is 
to drop the terms for which oj ^ -co', is equivalent to the Schrodinger picture 
interaction Hamiltonian: 

H I = Y 4»(w) ® B a (-u) = 5] Aa(u>) ® Bl(uj). (17) 

a,cu a,LU 

LEMMA 2. The Rotating Wave Approximation imply the conservation of the free 
energy of the global system, that is 

[H s + H B ,H]=0 (18) 
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2.1. Proof 

The necessary condition involved in the previous proposition is equivalent to 
the equation [Hs + Hb, Hi] = we are going to demonstrate. 

[H S + H B ,H] = [H S + H B ,H I } = [H S ,H I ] + [H B ,H I } (19) 

= -^2ujA a {oj) ®B a (-uj) +^2ujA a {uj) ®B a (-oj) = 0. 

a, a; a,co 

where we have made use of eq. (|9ll0p ■ 

LEMMA 3. The detailed balance condition in the thermodynamic limit imply U6\[ 

la p(uj) = e-^ la(3 (-uj) (20) 

where (3 = (ksT)^ 1 m 

COROLLARY 4. Let us suppose the temperature of the thermal reservoir to be 
the absolute zero, on the ground of Lemma 2 immediately we see that 

-f a p(u<0) = m (21) 

Let us now cast eq. ([T]) in a slightly different form splitting the sum over the 
frequency, appearing in eq. ([2]), in a sum over the positive frequencies and a sum 
over the negative ones so to obtain 

D( P s(t)) 

= £ 7 ^M(AV)p<^t M 

- l{Att(u,)A a (u>),p 8 }) 

+ ]T ^(-u)(A a Hu;)p s A^co) 

w>0,a,/3 

- ±{A*( U )Att( U ),p 8 }), (22) 

where we again make use of eq. (I13D . In the above expression we can recognize 
the first term as responsible of spontaneous and stimulated emission processes, 
while the second one takes into account stimulated absorption, as imposed by the 
lowering and raising properties of A a (to). Therefore if the reservoir is a thermal 
bath at T = the corollary 4 tell us that the correct dissipator of the Master 
Equation can be obtained by suppressing the stimulated absorption processes in 
eq. pa. 
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2.2. NuD Theorem 

We are now able to solve the markovian master equation when the reservoir 
is in a thermal equilibrium state characterized by T = 0. We will solve a Cauchy 
problem assuming the factorized initial condition to be an eigenoperator of the free 
energy H$ + H B - This hypothesis doesn't condition the generality of the found 
solution being able to extend itself to an arbitrary initial condition because of the 
linearity of the markovian master equation q 

NUD THEOREM 5. If eq. HJ) is the markovian master equation describing the 
dynamical evolution of a open quantum system S, coupled to an environment B, 
assumed to be in the detailed-balance thermal equilibrium state characterized by a 
temperature T=0, and if the global system is initially prepared in a state p(0) = 
PB (0)ps(0) so that (H s + H B )p(0){H s + H B ) = E 2 L p(0), where E L = E s + E B 
is the free energy of the global system then ps(t) is in the form of a Piecewise 
Deterministic Process J 13\j . that is a process obtained combining a deterministic 
time- evolution with a jump process. 

The proof of the theorem is contained in the paper [4j. My aim here is to give 
an explanation of the found implication. 

A PDP is a statistical mixture of alternative generalized trajectories evolving 
individually in a deterministic way. This statement is mathematically given by 
the equation 

N 

Ps(t)=Y,Pi( t )> ( 23 ) 

i=0 

where the quantum trajectories Pi{t) are obtained by the deterministic non- 
unitary equation 

Pl (t) = U(t)fi(t)tf(t), (24) 

where, in particular, 

f N (t) = PN (0) (25) 

and U(t) = e~i Bt , rf(t) = e^ Bh , B being 

B = H °-in £ 7^(«)^(«)A» = flb-^, (26) 



It is out of relevance to consider initial condition having non-zero coherence between the 
environment and the system because it is not possible to resolve them in the reduced dynamics 
obtained tracing on the environment degrees of freedom. 
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with H' hermitian. Finally, 



fN-j (t)= £ E - E 7 a ^)7^V)7^V0---7^V) (27) 

u>' ,a' ,f)' u" ,a" ,(3" wJ ,ai ,j3i 

x [[[...[ dt'dt" ...dtW^it^A 13 ' \u')U{t')U- x {e)A^' (w")tf (*")... 

JO JO JO JO 

u-\#)AP(ut)u{i?)f N {i?)rt{i?)A^\d)rt-\ir)... 

^ t (t")A Q "t( w »)^t-i(t'')t^(t')A a, t( £l /) ? 7t-i( t ') j j = i,...,at ■ 

These last are generalized respect to F.Petruccione and H.J.Carmochael ap- 
proach, which leads to ps(t) = Y^ilrfi >< V'il- The last expansion, in terms of 
proper trajectories, is obtainable from the mine if and only if we are able to put 
into diagonal form the spectral correlation tensor, that is known to be always pos- 
sible because of the positivity of 7°^, but nobody is able to do it, with exception 
of few highly symmetrical systems. 

The found solution (NuD theorem) ensures that the dynamical processes, 
whose statistical mixture gives the open system stochastic evolution, are deter- 
ministic. This demonstrates that the evolution is representable as a Piecewise 
Deterministic Process (PDP) [13]. The found solution generalizes the PDPs in- 
troduced by H.J.Carmichael and formalized by F.Petruccione and H.P.Breuer. 
Actually, it is applicable also when the Markovian Master Equation isn't in the 
Lindblad form. This, as already highlighted, in general, introduces simplification 
in the further calculations, but because of the difficulty to recast the equation in 
this form the results obtained are in general merely formal. Tough the eq. (|27p 
seems complicated to use it is a powerful predictive tool. We have tested it de- 
riving the photocounting formula [8], [20]; reproducing the environment-induced 
entanglement between two two-level not-direct-interacting atoms placed in fixed 
arbitrary points in the free space [U [9j [TOj \TT\ and Carmichael unravelling of the 
Master Equation [8| [T7]. 

Moreover I have tested the NuD theorem's predictive capability solving the dy- 
namics of two two-level dipole-dipole interacting atoms placed in fixed arbitrary 
points inside a single mode cavity in presence of atomic spontaneous emission and 
cavity losses [7]; n two- level not-direct-interacting atoms placed in fixed arbitrary 
points inside a single mode cavity in presence of atomic spontaneous emission and 
cavity losses [6]; a bipartite hybrid model, known as Jaynes-Cummings model, con- 
stituted by an atom and a single mode cavity linearly coupled and spontaneously 
emitting in the same environment (next subsection) and two harmonic oscillator 
linearly coupled and spontaneously emitting in the same environment (work in 
progress) . 

Two of my results are already published, others will be object of future papers. 
All of them contain the same (predictable? [4]) result: multipartite systems, 
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discarding the physical nature of the parts and of the environment, can exhibit 
entangled stationary states towards the system can be guided by a probabilistic 
scheme of measurement. 

3. Dissipative Jaynes and Cummings model 

The Jaynes-Cummings Model describes, under the Rotating Wave Approxi- 
mation, the resonant interaction between a single two-level atom and the single 
mode of the electromagnetic field protected by a perfect cavity (no loosing of en- 
ergy). The model has been introduced in 1963 [?] in order to analyze the classical 
aspects of spontaneous emission and to understand the effects of quantization on 
the atomic evolution. Actually, despite its apparent simplicity, this model has 
revealed interesting non-classical proprieties characterizing the matter-radiation 
interaction. Moreover, thanks to the recent experimental implementation of high 
Q cavities, it is today possible to verify the most of the theoretical predictions of 
the model [D [H H3] . 

The major experimental limitation is related to the coupling with a chaotic 
environment able to destroy the quantum coherences. A theoretical approach 
including the loss of energy due to the interaction of the atom and the cavity 
mode with the free electromagnetic field is more complete and, as we will show, it 
is suitable to reproduce the experimental measured decay of the population of the 
atom [lj. In particular, assuming a common bath of interaction between the cavity 
mode and the atom, the theoretical probability to find the atom in the excited state 
performs Rabi oscillation exponentially decaying. This fact is consistent with the 
open dynamics but it is not the only effect. Actually, the common bath induces 
cooperation between the two involved parts (mode and atom). This behavior 
competes with the exponential decay. In the long time limit the the exponential 
decay wins on cooperation if we work under the experimental condition performed 
by Haroche group. In the paper [1] is reported the experimental graph relative to 
the probability to find the atom in its excited state as a function of the time. We 
can interpret the upper part of the figure as the exponential decay and the lower 
part of it (increasing of probability) as the cooperation induced by the common 
reservoir. The new theoretical approach, here presented, is better than the usual 
one (two independent baths) because it is able to reproduce the experimental curve 
in a complete way. In fact the two bath approximation keeps account only for the 
dissipation meaning that the Rabi oscillation of the atomic population goes to 
zero every period characterizing the free dynamics of the bipartite system. In this 
case the cooperative part of the dynamics disappears: the two parts do not speak 
trough the common bath, the main behavior being the dissipation of energy in 
the reservoirs. Moreover it is possible to demonstrate that single bath approach is 
more general than the other including it as particular case. This fact is very well 
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understood if the parts are, for example, two or more atoms, in which case the 
cooperation is the maximum one if the distance among atoms is small enough and 
it reaches its minimum value when the distance goes to infinity [6] . In the last case 
the out diagonal terms of the spectral correlation tensor go to zero meaning that 
the parts see independent reservoirs. The lack of a microscopical derivation of the 
coupling constant of the mode with the electromagnetic field makes difficult the 
analytical derivation of an analogue relation in the case here analyzed. Despite 
this fact it will be shown that the single bath case is more general than the other 
one because the independent baths case does not reproduce the out diagonal terms 
giving a simplified Master Equation unable to reproduce part of the experimental 
measurements. 

The Hamiltonian describing the open system is 

H = h^S z + hujQa ] a + n(eaS + + e*a ] S-) + hY,Uj:b\{k)bx{k) (28) 

fc,A 
k,X k,X 

If we make the position 

H s = K^-Sz + hcuo^a + %{eaS + + e*Q f 5_), (29) 

H B = hJ2u; f :b{(k)b x (k), (30) 

k,X 



2 

Hj = ]T A a B a , (31) 



o=l 



where 



Sl =£bh 6 A$+0jSA 6t A®]» ( 32 ) 



fc,A 



B 2 = Y,{ s iM^ + %x bJ ^)] (33) 

fc,A 

A 1 = (S+ + S.), (34) 
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A 2 = (a f + a) (35) 

we can describe the reduced dynamics of the bipartite system at T = by a 
Master Equation of the standard form 

p s (t) = ~[H S ,ps(t)}+D(j, s (t)). (36) 

In this expression we have neglected the Lamb-Shift. This approximation is 
made possible because we have considered, ab initio, a direct linear static interac- 
tion among the parts respect to which the Lamb-Shift is negligible. In the above 
equation 



D(Ps(t)) = 1 1 > l (S^p s S + -±{S + S-,ps}) (37) 



+ 7 2 ' 2 (ap5a f - -{a<a,ps}) 



2 



+ l 1,2 {ap s S + - -{a ! S-,ps}) 



2 



2 ' 1 (S-psa f - -{S+a,p s }) 



' t"( , 



where 



T^^EK-aI 2 ^-^) (38) 



k.x 



l 2 ' 2 = ^J2\Hx\ 26 ^k-^ ( 39 ) 



A;. A 



7 1 ' 2 = * E R < s k9kx) 6 ^k - w) + ih* E Im ( s lx9kx) 6 ^k ~ w ) ( 4 °) 



fc,A fc,A 



7 2 ' 1 = 7T E jRe ( s L^A) <5 K ~ uj )~ ih7T E M*^**)^* - w ) ( 41 ) 

fc,A fc,A 

The master equation for p$ can be solved applying the NuD theorem to this 
case: 

n 

tt?(t) = !>(*)> ^ N ( 42 ) 

i=0 
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where n is the number of the excitation initially given to the system (ps(0) = 
Pn(0)) and i is the index giving the number of excitations characterizing every 
quantum trajectory. 

The trajectories evolve in time in accordance with 

Pi (t)=N(t)f t (t)N\t), (43) 

where fi(t) is given by eq. (|27p and N(t) = e~^Bt is a nonunitary temporal 
evolution operator, B being, in general, non-hermitian as it appears from the 
following equation: 

B = ~n S z + /i(Wa + MS+a + PcJS- - i— , (44) 

where 

(45) 

(46) 

(47) 

(48) 

Let us suppose the system in the initial state characterized by (n— 1) excitations 
in the cavity mode and the atom in its excited state | + >: 

/ 9 5 (0)=|(n-l),+ x(n-l),+ |, (49) 

then every quantum trajectory pi{t) belonging to statistical mixture characterizing 
the dynamical evolution of the system will have the form 

Piit) = Pi^C*) | (i - 1), + ><(*- 1),+ | +pi A (*) I (*-l),+ >< »,-|(50) 
+ P2,u(t) \i,-><(i- 1), + I +P2, 2l (i) \i, -><*,- | . 

The highest energy subspace (i = n) is easily solved and the block vector 
relative to this subspace has the form: 



n 


= OJQ 


-J7 1 ' 1 


n-~ 


= OJQ 


-* 7 2 ' 2 


AI 


= £ - 


-i 7 2>1 


P = 


e* - 


itf+y 



e -272,2(n-^)t-7M* 

f 

+ (| A n | 2 + | A | 2 ) cosh(26 n t) - 2ib n A(sin(2a n t) + sinh(26 n t))} 



P^(t) = : 7n~ri2 {(l^n| 2 -|A| 2 )cos(2a n t) (51) 



Pl,2 n (t) 



H^Ml e -2 72 , 2 (n-I)t- 7l , 1 t {sin(2ant) _ is[nH2bnt) (52) 



2A* n 



+ i——(cos(2a n t)—cosh.(2b n t))} 
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P2,i n (t) = P U n (t) (53) 



P2,2 n (t) = ?j f l | 2 e- 2 ^- 2 ("-^) t -^^{ S m(2a n t) -ismh(2b n t)}, (54) 
4 | A n | 



where 



A n = a n + ib n = \/A 2 + nMP (55) 

A = -(71,1-72,2)- (56) 

Let us note that if we have started from the initial condition | V'i(O) >=| n, — > 
we have obtained the same dynamical behavior. This fact ensures that an arbitrary 
linear combination of the two different initial condition will bring to the same 
dynamics. This fact is really important because it is not simple to prepare one or 
the other of the initial states. Actually, when we inject an excitation inside the 
system we can only know that the system is in a statistical mixture of the two 
states. But we have seen seen that the dynamics is not case sensitive and therefore 
a statistical mixture of the two states leads to the described dynamical evolution. 

3.1. Entanglement building up 

The circumstance that we succeed in finding the explicit time dependence of 
the solution of the master equation (|36|) provides an occasion to analyze in detail 
at least some aspects of how entanglement is getting established in our bipartite 
system. As particular case we can choose n = 1 so obtaining in a simple way the 
complete dynamics of the open system in the form 

Ps(t) = P o(t) + pi(t), (57) 

where 

Pi(t) = pi,i 1 (t)|0,+ ><0,+ |+pi A (t)|0,+ ><l,-| (58) 

+ /02,1a (t) I 1, - >< 0, + I +p 2 , 2l (t) I 1, - >< 1, - I 

and 

po(t) = (1 - pi M (t) - p 2 , 2l (*)) I 0, - X 0, - I (59) 



On the basis of the block diagonal form exhibited by Eq. (|57p , at a generic time 
instant t, the system is in a statistical mixture of the vacuum state of the system 
and of a one-excitation appropriate density matrix describing with certainty the 
storage of the initial energy. In order to analyze the time evolution of the degree of 
entanglement that gets established between the two initially uncorrelated parties, 
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Fig. 1: Blablabla... 



we exploit the concept of concurrence C first introduced by Wootters [HI S3] . If, 
at an assigned time t, no photon have been emitted the conditional concurrence 
C assumes the form: 



C{t) 



/ 2(pi,l n p 2 ,2 n + | Pl,2 n I 2 ) +4/?2,2„ | Pl,2 n 



Pl,l„ + P2,2 n 



(60) 



/ 2(pi,l n p 2 ,2 n + | Pl,2 n | 2 ) -4/) 2 ,2 n | Pl,2„ 



V Pi, In + P2,2 n 

In the analyzed case (n = 1), as clearly shown in Fig.l, obtained using the ex- 
perimental values setted by Haroche's group, the degree of entanglement starting 
from zero increases during the transient collapsing to the initial value when time 
is long enough. This fact depend on the choice of the atom whose spontaneous 
emission time 711 is much longer than the cavity damping time 72.2- In accor- 
dance to this fact the probability to find the atom in the excited state starting 
from 1 go to zero when t 3> (71 1 + 72,2)"* as clearly showed in Fig. 2. This Figure 
reproduce in a perfect way the experimental measures performed by Haroche's 
group [1J. The standard theoretical models assume two different channel of dissi- 
pation (one for the atom and one for the cavity mode). The corresponding Master 
Equation is simpler to solve because of the absence of cooperative terms |3j but the 
corresponding dynamics fits only the upper part of the measures of the Haroche's 
group {dissipative behavior). The low part of the graph represent the cooperation 
induced by the common reservoir between the cavity mode and the atom. 

Such cooperation become the maximum one when 7^1 = 72,2 and fase, as 
clearly shown in Fig. 3 and Fig. 4. These ones depicted the concurrence and the 
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probability to find the atom in its excited state, respectively. 

Under this condition (Decoherence Free Regime) Eq. (|57p suggests that, for 
t S> (71^ + 72,2) ~ 1 , the correspondent asymptotic form assumed by ps(t) is time 
independent and such that the probability of finding energy in the bipartite system 



is 



2" 



ps(t> (71,1 +72,2) ) 



+ - 



1 

2 

1 



ipG >< ^G 
IpA >< tpA 



(61) 



where 



I ipG >=\ 0,- >< 0, 
is the ground state of the bipartite system and 

1 



^4>=-y=(|0,+ Xl, 



1,- ><0,+ 



(62) 



(63) 



is the maximally antisymmetric entangled state of the system. 

This fact suggests that stationary entangled states of the JC system can be 
generated by putting a single photon detector able to capture in a continuous way 
all the excitations lost by the system in the reservoir. Reading out the detector 
states at t 3> (71,1 + 72, 2) -1 ; if no photons have been emitted, then, as a conse- 
quence of the measurement outcome, our system is projected into the maximally 
antisymmetric entangled state | ip >= -4=(| 0,+ >< 1, — | — | 1,— >< 0,+ |). 
This is the main result of the paper which means that a successful measurement, 
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performed at large enough time instants, generates un uncorrelated state of the 
two subsystems, bipartite system and reservoir, leaving atom and cavity in their 
maximally antisymmetric entangled state. 

This ideal result has to be corrected by the introduction in the microscopical 
model of a second bath of interaction able to take account of the cavity leakage 
of energy because of the imperfect mirrors. In terms of the hamiltonian operator 
this means: 

H = n'^-S z + huJoa t a + h(eaS + + e*a^S-) + fiY / ^kb[(k)b x (k) (64) 



2 



k,\ 



+ fi£«jcl(fe> A (fel + +5>| x b x (k) + g lx b\(k)](S + + S.) 



k,X k,X 



k,X k,X 

If we make the position 



Hs = h — S z + fiuoora + h(eaS + + e*cv S-), (65) 



H B = hJ2^b{(k)b x (k) +hJ2^c{(k)c x (k), (66) 

k,X k,X 

2 

H I = Y / A a B a , (67) 

01=1 



where 



Bl = Hi9lM^ + 9k,x b ^(k)}, (68) 
k,x 

b2 = Y} S \M^ + %x b ^(k)] + J2i r lx c ^) + r^c\{k)] (69) 

k,X k,X 

A 1 = (S+ + S.), (70) 

A 2 = (^ + a) (71) 
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we can describe the reduced dynamics of the bipartite system at T = by a 
standard Master Equation and we can solve it in the same way of the previous case. 
The changes in the microscopical model do not introduce variation in the formal 
solution. Instead, the presence of two different channel of dissipation modifies the 
dynamical behavior. Actually, the system has now the possibility to loose energy 
in environments that do not speak each other. This means that, when the time is 
much longer than the sum of the single emission time, the coherence induced by the 
common bath during the transient will go to zero in the long time domain. Despite 
this fact, the dechoerence time can be made as long as we need to implement the 
required quantum protocol. Actually, named k the cavity decay rate, if k is much 
greater than 7^1 = 72,2, then the storage of energy can be maximized for a time 
sufficient to realize the quantum protocol. 

4. Conclusion 

In this paper we have considered the interaction of a Jaynes and Cummings 
system with the electromagnetic field (and with another phenomenological zero 
temperature bath) in its vacuum state and, solving the dynamical problem, we 
have analyzed the amount of entanglement induced in the bipartite system by the 
common electromagnetic reservoir. This has allowed us to quantitatively charac- 
terize the regime under which field-induced cooperative effects are not vanished by 
dissipation. Once the Decoherence Free Regime is reached, transient entanglement 
tends to become stationary and, therefore, usable for quantum gate implementa- 
tion. 

The asymptotic solution of the dynamical problem appears to be a statistical 
mixture of a maximally entangled state and the ground state of the open system, 
the probability to obtain one or the others being the same. In the whole temporal 
domain the found solution tell us that the state of the system is a statistical mix- 
ture of the free energy system eigenoperators. This fact is general enough [4] and it 
is consistent with the existence of a photon detector device because the act of mea- 
surement introduces a stochastic variable respect to which we can only predict the 
probability to have one or another of the possible alternative measures [4] . These 
probabilities can be regarded as the weight of the possible alternative generalized 
trajectories. With this approach the dynamics has to be depicted as a statistical 
mixture of this alternative generalized trajectories. Moreover the found trajecto- 
ries evolve in time in a deterministic way: for example the trajectory relative to the 
initially excited system state is a shifted free evolution characterized by complex 
frequencies that means an exponential decay free evolution. This statement may 
give the sensation that every system has to decay in its ground state because of the 
observed dynamics. It is in general not true. Actually, if the system is multipartite 
as ours, it is possible that it admits excited and entangled equilibrium Decoher- 
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ence Free Subspace (DFS) [l6](so as it happens in some highly symmetric models), 
constituted by states on which the action of Hj is identically zero. If the system, 
during evolution, passes through one of these states, the successive dynamics will 
be decoupled from the environment evolution. An equilibrium condition is reached 
in which entanglement is embedded in the system. 

Bibliography 

1. J.M.Raimond, M.Brune, S.Haroche, Rev. Mod. Phys., 73, 565 (2001) 

2. T. Nirrengarten, A. Qarry, C. Roux, A. Emmert, G. Nogues, M. Brune, J.-M. Raimond, and 
S. Haroche, Phys. Rev. Lett., 97, 200405 (2006) 

3. G.M. Palma, K-A Suominen, A. Ekcrt, Quantum computers and dissipation, Proc. Roy. Soc. 
Lond. A 452 567-584 (1996) 

4. S.Nicolosi, Open Systems and Information Dynamics 12, (2005) 

5. E. Knill, R. Laflamme, H. Barnum, D. Dalvit, J. Dziarmaga, J. Gubernatis, L. Gurvits, G. 
Ortiz, L. Viola, W.H. Zurek \ quant-ph/0207171\ 



6. S. Nicolosi, A. Napoli, A. Messina, quant-ph0306104 (2003), submitted to EPJD 

7. S. Nicolosi, A. Napoli, A. Messina , F.Petruccione quant-ph0402211 (2004), axcepted on 
PRA 

8. S.Nicolosi, Soluzione operatoriale esatta di master equation ottiche generalizzate: appli- 
cazioni a problemi di interazione radiazione materia, Thesis (2003), unpublished 

9. Plenio, Huelga, Beige, Knight, Phys. Rev. A 59 (1999) 2468 

10. T. Pellizzari, S.A. Gardiner, J.I. Cirac, P.Zoller, Phys. Rev. Lett. 75, 3788 (1995) 



11. Z.Ficek, R.Tanas, quant-pi 0302124 

12. A.Olaya-Castro, Niel F.Johnson, quant-ph0306104 (2004), to appear in Handbook of Theo- 
retical and Computational Nanotecnology 

13. F.Petruccione, H.P.Breuer, The Theory of open quantum sistem, Oxford University Press 
(2002) 

14. C.W.Gardiner, Quantum Noise, Springer- Verlag Berlin (1991) 

15. W.H.Louiselle, Quantum statistical properties of radiation, Ed. John Wiley e Sons (1973) 

16. R.Alicki, K.Lendi, Quantum Dynamical Semigroups and Applications, Lectures Notes in 
Physics V286, (Springer- Verlag 1987) 

17. H.Carmichael, An Open System Approach to Quantum Optics, Lectures Notes in Physics 
ml8, (Springer- Verlag 1993) 

18. Y.Castin, K.Molmer, Phys. Rev. Lett. 74, 3772 (1995) 

19. J.Dalibard, Y.Castin, K.Molmer, Phys. Rev. Lett. 68, 580 (1992) 

20. E.B.Davies, Commun. Math. Phys. 15, 277 (1969) 

21. M.D.Srinivas, E.B.Davies, Opt. Acta 28, 981 (1981) 

22. M.B.Plenio, P.L.Knight, Rev. Mod. Phys. 70, 101 (1998) 

23. F.Petruccione, H.P.Breuer, Concepts and methods in the theory of open quantum systems 
systems, Lectures Notes in Physics V622, Springer- Verlag Berlin (2003) 

24. H.J. Carmichael and K. Kim, Optics Commun. 179, 417 (2000) 



[Author and title] 20 

25. J. P. Clemens, L. Horvath, B.C. Sanders, and H.J. Carmichael, Phys. Rev. A 68, 023809 
(2003) 

26. V.B.Braginsky, F.Ya.Khalili, Quantum Measurement, Cambridge University Press (1992) 

27. S.M.Barnett, P.M.Radmore, Methods in Theoretical Quantum Optics, Oxford Series in Op- 
tical and Imaging Science, 15 (2002) 

28. C.Leonardi, F.Persico, G.Vetri, La rivista del Nuovo Cimento9(1986) 3,4 

29. G.Benivegna A.messina, J.Mod.Opt. 36 (1989) 1205 

30. M.Kozierowski, S.M.Chumakov, A.A.Mamedov, Journal of Mod. Opt. 40 (1993)453 

31. B.W. Shore, P.L. Knight, Journal of Mod. Opt. 40 (1993)7 

32. I.K.Kudryavstev, A.Lambrecht, H.Moya-Cessa, P.L.Knight, Journal of Mod. Opt. 40 (1993)8 

33. S. J. D. Phoenix, S.M.Barnett, Journal of Mod. Opt. 40 (1993)6 

34. Dicke, phys. rep. 93 (1954)99 

35. M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cam- 
bridge University Press, Cambridge, (2000). 

36. Guo Ping Guo et al, Phys. Rev. A 65, 042102 (2002) 

37. Lu-Ming Duan et al, Phys. Rev. A 58, (1998) 

38. J.I.Cirac, P.Zoller, Phys. Rev. A 50, R2799 (1993) 

39. Shneider, Milburn, Phys. Rev. A 65, 042107 (2002) 

40. G.J.Yang, O.Zobay, P.Meystre, Phys. Rev. A 59, 4012 (1998) 

41. E.Hagley, X.Maitre, G.Nogues, C.Wunderlich, M.Brune, J.M.Raimond, S.Haroche, Phys. 
Rev. Lett. 79, 1 (1997) 

42. P.Foldi, M.G.Benedict, A.Czirjak, Phys. Rev. A 65, 021802 (2002) 

43. M. Brune, F. Schmidt-Kaler, A. Maali, J. Dreyer, E. Hagley, J. M. Raimond, S. Haroche 
Phys. Rev. Lett. 76, 1800 (1996) 

44. S. Hill, W. K. Wootters Phys. Rev. Lett. 78, 5022-5025 (1997) 

45. W. K. Wootters Phys. Rev. Lett. 80, 2245-2248 (1998) 



46. D.A.Lidar, B.Whaley quant-ph/0301032 



